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We present calculations of matrix elements of 4-quark operators in the pion and in the nucleon extracted from 
quenched Monte Carlo simulations at /3 = 6.0 using Wilson fermions. These operators are relevant for higher-twist 
effects. We are particularly careful to avoid mixing with lower-dimensional operators by choosing appropriate 
flavour structures. 



1. INTRODUCTION 

The operator product expansion (OPE) ex- 
presses (Nachtmann) moments of hadronic struc- 
ture functions in terms of Wilson coefficients 
(usually calculated in perturbation theory) and 
nonperturbative hadronic matrix elements. The 
leading contribution in the deep-inelastic limit 
Q 2 — y oo is provided by operators of twist 2 with 
corrections suppressed by powers of 1/Q 2 coming 
from operators of twist 4 and higher. Schemati- 
cally one finds for n = 2, 4, 6, . . . 



is t — 2. In the flavour-nonsinglet channel, the 
twist-2 operators are 2-quark operators, 

^ 1 D fl2 ---D^, (2) 

symmetrised in all indices and with trace terms 
subtracted. 

An important class of twist-4 operators are 4- 
quark operators. In particular, the twist-4, spin- 
2 matrix element A^' is given by 

(p\A c { ^ } - traces |p) = 2A ( ^ ) (p f _ l p l/ - traces) (3) 
in terms of the 4-quark operator 



dxx n - 2 F 2 (x,Q 2 ) 



Nachtmann 



(4) 



c£ 2) (q 2 /m 2 ,s(m))4 2) M 



(i) 



, ci 4) (Q 2 /M 2 ,fl(M)) , W 
Q 2 



The reduced matrix elements A$ of twist t and 
spin n depend on the renormalisation scale //. 
Whereas the corresponding Wilson coefficients 
c£} are dimensionless, the mass dimension of An 
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The quark field ip carries flavour, colour, and 
Dirac indices, the matrices t a are the usual gen- 
erators of colour SU(3), and the flavour matrix G 
contains the quark charges: 

G = diag(e u , e d ) = diag(2/3, -1/3) (5) 
for two flavours. The Wilson coefficient reads (!]] 

4 4) = 3 2 (i + o( 3 2 )). 

These expressions are to be compared with 
their twist-2 counterparts: 

(p\0{nu} - traccs|p) = 2A\ '{p^Pu - traces) (6) 
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with the operator 
and the Wilson coefficient c 



(7) 



(2) 



l + 0( 5 2 ). 



The operators (Q) and (0) transform identically 
under Lorentz transformations, but (^) has di- 
mension 6, wheras (^) has only dimension 4: 4- 
quark operators will in general mix with 2-quark 
operators of lower dimension. This fact compli- 
cates the investigation of 4-quark operators, be- 
cause the mixing with lower-dimensional opera- 
tors cannot be calculated reliably within pertur- 
bation theory. For the time being, we do not at- 
tempt a nonperturbative calculation of the renor- 
malisation and mixing coefficients of 4-quark op- 
erators. Instead we restrict ourselves to cases 
where mixing with lower-dimensional operators 
is prohibited by flavour symmetry. 

In the following we present Monte Carlo data 
from quenched simulations at f3 = 6.0 with Wil- 
son fermions on a 16 3 x 32 lattice. From now on 
all operators are written down in Euclidean space. 

2. PION 

In the case of the pion, we consider the sym- 
metry SU(2)p, i.e. isospin symmetry. While 2- 
quark operators can have at most isospin I = 
1, there are 4-quark operators with I = 2. 
Such operators are protected from mixing with 
lower-dimensional operators and it makes sense 
to renormalise them perturbatively. Therefore we 
study the 1 = 2 pion structure function 



Ft 2 = Ff 



-2F? 



(8) 



which is purely higher twist and receives no con- 
tributions from 2-quark operators. 

Omitting Dirac and colour matrices, the flavour 
structure of the relevant operator in the OPE is 
given by 

(e u uu + eddd)(e u uu + eddd) 



-\{uu)(uu) + (dd)(dd) 
6 



(uu)(dd) — (dd)(uu) — (ud)(du) — (du)(ud)] 



up to contributions with / = 0, 1. Although this 
flavour structure excludes mixing with 2-quark 
operators, mixing with other 4-quark operators 
is still possible. For spin 2, a basis of 4-quark op- 
erators whose 1 = 2 components are closed under 
renormalisation is given by the spin-2 projections 
of 



va 



(9) 



together with the operators V^„, A^, and T^ v 
which differ from the above by the omission of 
the colour matrices t a . One-loop results for 
the renormalisation coefficients are given in [Q]. 
These are to be combined with the bare matrix el- 
ements extrapolated to the chiral limit. Express- 
ing the matrix elements in terms of the pion decay 
constant f n we get H 



.4. 



(4) 7=2 



0.133(51) / T 2 , 



(10) 



where = 131MeV in the real world. The order 
of magnitude of this result is reproduced by the 
vacuum insertion approximation, which yields an 
answer proportional to f 2 . For the first moment 
of F2 we find 



f 1 

/ dxF 2 (x,Q 2 )\^ 
Jo 



Nachtmann 



(ii) 



1.67(64) ■ 



Q 2 



Q(a 2 s 



Except for very small values of Q 2 this is consid- 
erably smaller than the leading twist-2 contribu- 
tion to the 7r + structure function computed in a 
quenched Monte Carlo simulation ||: 



dxF 2 (x,Q 2 )\l 



Nachtmann 



0.152(7). 



(12) 



Of course, the fact that the particular twist-4 
contribution ( |li"| ) is rather small does not tell us 
that all other higher-twist contributions will be 
small as well. 
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3. PROTON 

In the case of the proton it is more difficult to 
find 4-quark operators that are safe from mixing 
with lower-dimensional operators: the expecta- 
tion value of any 1 — 2 operator in the proton 
vanishes. So we have to enlarge the flavour sym- 
metry which we consider from SU(2)p to SU(3)f, 
i.e. we assume three quarks of the same mass. 
Correspondingly, the flavour structure of the op- 
erator in the OPE is now 

O = (e u uu + eddd + e s ss){e u uu + eddd + e s ss) . 



Whereas 2-quark operators 
SU(3)f according to 3 ® 3 = 
for 4-quark operators: (3 C8 
2-l©4-8ffil0ffil0ffi27. 



transform under 
= 1 © 8, we have 

3) <g> (3 <8> 3) = 
4-quark operators 



with J = 0, 1, J3 = 0, and hypercharge Y = 
from the 10, 10, 27 multiplets do not mix with 
2-quark operators and can be used in a proton 
expectation value, e.g. the 1=1 operator 



o I=1 



1 



27 = Y^[i u u)(uu) - (dd)(dd) 



— (us)(su) — (su)(us) + (ds)(sd) + (sd)(ds) 

— (ss)(uu) — (uu)(ss) + (ss)(dd) + (dd)(ss)] 

belongs to the 27 multiplet. Analogously to the 
case of the pion, one can find linear combinations 
of the structure functions of the octet baryons (p, 
n, A, E, S) which project out the desired flavour 
component, e.g. 

(S+|0|E+) - 2(S°|0|S } + (E-|0|E-> 



(£+|0|£+) + (E-|C?|£-) 



6(A|G|A) - 2(H U |0|H°) - 2<S-|0|S-) 



2{p\0\p)-2(n\0\n) = -(p\0 J 2 ^\p) . 



The proton matrix elements are computed 
in the standard fashion from ratios of 3-point 
functions (B(t)O(T)B(0)) over 2-point functions 
(B(t)B(0)) (0 < t < t) with suitable interpolat- 
ing fields B and B. For a general 4-quark opera- 
tor the 3-point function (B(t)O(T)B(0)) consists 
of three types of contributions, which can be rep- 
resented pictorially by the following diagrams 



OO 





It is precisely through contributions of the form 
of the first two diagrams that the mixing with 
lower-dimensional operators occurs. Therefore 
these contributions cancel in the operators which 
we consider, and we are left with the contribu- 
tions of the last type only. 

As an example of our results we show in fig]l] 
the chiral extrapolation of the bare proton ma- 
trix element of V£ 4 — traces (1=1 component in 
the 27 representation of SU(3)f) divided by the 
fourth power of the proton mass m p . It indicates 
the order of magnitude of all 4-quark matrix el- 
ements that we have studied. After renormali- 
sation this expression contributes to the reduced 
matrix element A^f 1 according to 



.4 



(4) 



in- 
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(13) 



and we obtain for the lowest moment of F% in our 
special flavour channel 



/ dxF 2 {x,Q 2 
Jo 



1 27,7=1 

I Nachtmann 



-0.0006(5) p ' + 0{a 2 s ) 



(14) 



In the proton the corresponding twist— 2 contribu- 
tion is w 0.14. Once again, the twist-4 correction 
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Figure 1. Chiral extrapolation of the bare proton 
matrix element of O^f 1 for the operator U 4 C 4 — 
traces, divided by m^. 



is tiny. This may be compared with bag model 
estimates, which give numbers oc B/rrip « 0.0006 
for the prefactor in dlj), where B w (145MeV) 4 
is the bag constant 

4. COMPARISON OF RESULTS 

Let us finally compare in fig]^ the renormalised 
pion matrix elements (7r + |6 • (• • ■) i=2 \tt + ) /m 2 
with the corresponding renormalised proton ma- 
trix elements (p|10 • (• • •)27 =1 |p)/ m p ( m lattice 
units). In the operators we have set fi = v = 4 
(with traces subtracted). The normalisation of 
the pion and proton matrix elements is chosen 
such that the flavour structure (uu)(uu) appears 
with the factor 1 in both cases. In view of the 
fact that pion and proton are very different par- 
ticles it is hardly surprising that the numbers do 
not show many similarities. 
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Figure 2. Renormalised 4-quark matrix elements 
in the pion and in the proton. 

Unfortunately, these flavour structures are some- 
what exotic and not easily accessible to experi- 
ment. The investigation of physically more inter- 
esting flavour channels is of course desirable, but 
impossible without progress in nonperturbative 
renormalisation. It remains to be seen what will 
happen first: a measurement of structure func- 
tions in the above flavour channels or a reliable 
lattice calculation of twist-4 matrix elements rel- 
evant for the proton or the pion. 
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5. SUMMARY 

Contributions of higher twist to hadronic struc- 
ture functions promise to challenge lattice QCD 
for a few more years. Still, 4-quark operators 
can give reasonable signals in present quenched 
Monte Carlo simulations. However, when the 
flavour structure is such that it prohibits mix- 
ing with lower-dimensional (2-quark) operators, 
the matrix elements turn out to be rather small. 
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